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1 Introduction 



The notion of light-cone wave functions and distribution amplitudes is a well 
known tool to study the physics of light mesons[|l], ||, [§. These concepts were 
indroduced in the framework of light-cone quantification; they also appear 
in the parametrization of non-perturbative behavior of light mesons in the 
infinite momentum frame [3J. 

The different distributions amplitudes are the basic non-perturbative ob- 



jects needed for the QCD light-cone sum rules ||10|| , which is one of the meth- 
ods to compute exclusive decays of heavy mesons into light mesons (an im- 
portant way to test the Standard Model and its extensions). In order to 
obtain more precise results, the distribution amplitudes generate by a non- 
zero meson mass is usefull. 

The most popular model for these distribution amplitudes is the con- 
formal expansion, where they are developped on a (almost) multiplicatively 
renormalizable basis [[|, ||. In this framework, the meson mass corrections are 
obtained via the exact equation of motion!!. But other models exist (renor- 



malon model in |DJ, instanton model in therefore a model-independant 



way to compute the meson mass corrections is needed. 

The method used in this letter consist of expressing the distribution am- 
plitudes generated by the mass of a vector meson in term of the leading twist 
distribution amplitude, keeping only the (next-to-) leading terms in the Op- 
erators Products Expansion on the light- cone 0]. This idea was introduced 
in || in the case of a scalar theory. 

The paper is organized as follows. After the introduction of the usefull 
notations, the different notions of twist is explained, on which the method of 
computing mass correction is based. 

Then the case of chiral even and chiral odd two-points distribution are 
treated. 

2 Notations 

The two-points distrubtion amplitudes considered are extracted from this 
kind of matrix element : 

(0\u(x)[x, -x)Td(-x)\p) (2.1) 
Where x is almost on the light cone, u G [0,1], T any kind of product 
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of 7 M matrices and [x, y] a path-ordered gauge factor along the straight line 
connecting x and y : 



[x, y] = P exp 



dt{x - y)^{tx + (1 - t)y) 



(2.2) 



This gauge factor is a way to introduce interaction (see [JP2]]), it ensure 
the gauge invariance of these non-local matrix elements (I will omit to write 
it sometimes). 

These matrix elements depend on three vectors : 

- Pfj, : momentum of the p-meson 

- : polarization vector of the p-meson 



x. 



With the relations : 



p2 

e (A) 

P • c w 



e (A) . e (A) 



-1 




(2.3) 
(2.4) 
(2.5) 



3 Definitions of twist 

The different contributions to the matrix elements |2.1| are classified by their 
"twist", but there are two different definitions of this concept: 

For local operators, the "twist" means dimension minus Lorentz-spin. 
To extract the contribution of a definite twist for non-local matrix elements 



like |2.1| , an expansion in local operator (OPE) has to be done and all local 
operators which have a definite twist in the sense given above must be re- 
summed. This will give the non-local term which has a definite power in x 2 . 
To understand how this definition of "twist" means "power in x 2 ", see J/J. 
From now on, this definition will be called "theoretical twist". 

The second definition is built in analogy to partons distributions. In the 
latter case, the different structures of non-local matrix elements are separated 
with their different powers of Q (the hard momentum transfer) in deep in- 
elastic scattering : a term of twist t contributes to the inclusive cross section 
with coefficients which contain t — 2 or more power of 1/Q. A good descrip- 
tion of this classification can be found in ||. From now on, this definition 
will be called "physical twist". 
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In this letter, the distributions amplitudes which are generated by meson 
mass corrections (these are non-leading "physical" twist distribution ampli- 
tudes) are computed in the leading (and next-to-leading) "theoretical" twist 
expansion. 



4 Meson mass correction 
4.1 chiral-even case 

The mass of the p-meson induces some non-leading (physical-) twist structure 
for the non-local matrix element: 



(0\u(x)^[x, -x]d(-x)\p-(P,\)) = 

■ x 



= f P m p 



+ 



M P-x 

e (A) • x 



1 / 2 2 

J due* p - x ( B ^{u) + X -^B^\u) 



1 e (A) -x . rl 
2 X ^~[P~xf m ' p 



due^ p - x C(uY e) 



(4.1) 



with £ = 2u — 1. 

Some wave functions would disappear if the mass of the p were 0. Hence 
mass corrections are higher (physical-) twist contributions to non-local matrix 
elements (A^ e \u) and C^ e \u) are twist 4, B^ is twist 5). 

The idea is to compute the mass corrections in the leading theoretical 
twist approximation. In this part, the case with the two-point chiral even 
distributions is computed. The "brute force" method, i. e. the expansion of 
the distributions in moments (x — > 0) can be used. It gives : 



(0\u llx (iDx) n d\p-(P,\)) 
e (A) • x 



= P 



" P-x 



(Px) n M^ 
i ■, , e (A) • x 



2 2 
X TYl 

—JL n {n-\)Mt^) n ' 2 



2 2 
X 771 

(Px) n M* - —^n(n - l)M*_ 2 (Px) n - 2 



1 e 
x 



(A) . 



X 



2 ^(P- 



x 



2 m 2 p (PxrM c n 



(4.2) 
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where the moments M£, M^, and Mjf are defined like in the fol- 

lowing way: M n = f Q du^ n f(u). 

The leading (theoretical-) twist of the operators in the left hand side can 
be written. The condition of symmetry and zero trace gives : 



(0\u^(iDx) n d\p-(P,\)) = 

1 " p-x [ } + n+ 1 v " p-x y 



(Pl) "kw " f'M^wl <<0 " >> 



The comparison of these two equations implies : 

Mi = <(O n » (4.4) 

< = -^-r«O n )) (4.5) 

71 + 1 

M t, = ( 4 - 6 ) 

< = (4-8) 

There is another method to get the leading (theoretical-)twist part of 
the distributions, giving directly integrated equations between the different 
wave functions. The different relations for non-local operators in appendix |A| 
are used: they give integrated forms for the leading twist part of two-point 
operators. 

One takes equation |A.4 : 

r 1 d _ 

P'(x)7 a ^(-x)] 8ym = J du—ip(ux)xip(-ux) (4.9) 

where x = x Ur y u . The subscript "sym." means that in the Operator Product 
Expansion (local expansion around x — 0), only the operators which have 
symmetrical Lorentz- indices are taken. But for the leading theoretical twist, 



(4.3) 
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the traces has also to be removed, so the leading twist part of the right hand 



side of 4.9 is also taken : 



(4.10) 



The condition for the leading twist part of ip{ux)xip{—ux) is the following 
(equation |A.5|) : 

'> " ' (4.11) 



Now the matrix element of ip(ux)xip{—ux) can be taken. The parame- 
terization of the leading twist part up to m 2 corrections is the following: 



(0| — x]d{— t \p (P,X)} 



(e (A) • x) I dugiu) 
'o 



1 1 

x m. 



+ —^f^p x ) + 0( 



(4.12) 



where g(u) is the leading twist wave functions in the limit m p — > 0. 



The condition [4. 1 1| can be applied; it gives (after some calculation) : 

f dug{u) [3f^Px) + (Px)f^Px) - i 2 e liPx ] = (4.13) 
Jo 

The following ansatz for f^(Px) can be done (having in mind the case of 
scalar field |§): 



k = a(0 / dvb(v) 



Then equation 4.13 gives 



b(y) = v 2 
This can be written in this form : 

{Q\[u{x)x[x, —x\d{~x)] l t \p~(P,\)) 
dug{u) [{e^-x)e* p \ 



(4.14) 



(4.15) 
(4.16) 



(4.17) 
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with 



[(eW-x)e^] L , 



3 (A) 



x) 



The equation [4 . 1 0| can be applied. After some algebra, one gets the 
distributions A^(u), B^(u), C (e) (u), i (e) (w) and B^(u) defined in |0 
the leading (theoretical-) twist approximation : 



in 



duB ie \u)e^ Px 
duC {e \u)e iiPx 

duM e) (u)e^ Px 



1. t. 



l. t. 



l. t. 



l. t. 



l ,i 

dug(u) / dte^ Px 
Jo 



(4.19) 
(4.20) 



dug{u) 



e^ Px -3 I dte^ Px + 2 / dt \ dve ivt ^ Px 

o </o 



dug{u)^ 2 I dtt 



dug{u) 



e itSPx_ 2 / dvv 2jvtdPx 



1 dti 2 f 1 v 2 dve** Px 
Jo 



(4.21) 
(4.22) 

(4.23) 



There is no index "1. t." for A^ e '(u) because this distribution contains only 
the leading twist part of the matrix element (0\u(x)'j f j i [x, — x]d(— x)\p~(P, A)). 
So these last equations show that all the distributions can be computed 
when one has A^{u), in the leading (theoretical-)twist approximation. They 
reproduced exactly the relations between moments (equations O to 
The advantage is that they are directly in an integrated form. 



4.1.1 Adding twist 3 and twist 4 terms 

The idea is to add contributions of (theoretical-)twist 3 and 4 without the 
gluonic contribution. In the appendix [A], there is the following operator 
relation ( A.lOj) : 

- - 1 f 1 - 

ip(x)xip(—x) = [tp(x)xip(—x)'j l — -x 2 / dvv 2 d 2 ip(vx)xtjj(—vx) 

4 Jo 

+ operators with gluons + . . . (4.24) 
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The matrix element of the term with a total derivative, gives 



x I dvv (0|<9 u(vx)x[vx, — vx]d(— vx)\p(P, A)) 



m 2 -x 2 I dvv 2 (0\u(vx)x[vx, — vx]d(— vx)\p(P, A)) (4.25) 



o 

1 

The 0(m°) term (equation |4.12| ) in the parameterization of the matrix 
element in the right hand side implies : 

I f 1 

— -x 2 / dvv 2 (0\d 2 u(vx)x[vx, — vx]d(— vx)\p(P, A)) = 

4 In 



;r) 



2 2 »1 



/ dug{u) [ dvv 2 e lv ^ Px (4.26) 
Jo Jo 



So the leading twist part with the total derivative term is : 
(0| [u{x)x[x, -x]d(-x)] L t + tot . dcr . \p~{P, A)) = 

dug(u) [(e^.x)e* p % t 



tot. der. 



(4.27) 



with 



[(eW - x)e* p *] L t 



tot. der. 



fe (A) • x) 



2 2 



(£ 2 + l) / dvv 2 e iv Z Px 
Jo 



(4.28) 



This last form can be used to give some twist 4 corrections to equation 



4.19 to 4.23 



Now the twist 3 part without gluon can be added. The matrix element 
of the equations |A.7| is taken, suppressing the gluonic terms : 



(0| [u{x)^[x, -x]d{-x)] t3 \p (P, A)) 



/ dttx u (-iP a )(0\u(txhf3^[tx,-tx}d(-tx)\p-(P,X)) 



(4.29) 



[u{x)j^ 5 [x, -x}d(-x)] t3 \p (P, A)) = 

/ dttx u {-iP a )(0\u{tx)j l3 [tx, -tx}d{-tx)\p~ (P, A)) 



(4.30) 



Combining these two equations, replacing the right hand side of the sec- 
ond one with the parameterization of |4.1| , one gets: 

(0| [u(x)^[x, -x]d(-x)] t3 \p~(P, A)) = 



(e^(P -x) 2 +P fl (e^ ■ x){P ■ x) f du [ dtt [ dvv 2 B ( - e 

Jo Jo Jo 

+ ml I du I dtt [ dvv 2 e ivt t px x 



e \u)e ivt ^ Px 



JO 



(e^(P ■ x) 2 + P M ( e W . X )(P . X )t^-B^ (u) + (x 2 e^ — x^e^ ■ x)) B^ e \u) 

(4.31) 

Including the leading twist term, the total derivative term ( |4.28| ) and this 
last equation, equations for A^ e \u), B^ e \u), C^ e \u), A^ e '(u) and B^ e \u) can 
be obtained: 

A^(u) = g(u) (4.32) 

/ duB (e \u)e iiPx = [ dug{u) [ dt^ Px 
Jo Jo Jo 

-{P-xf I du f dtt [ dvv 2 e mt ^ Px B (e) 
Jo Jo Jo 



U) 

(4.33) 



1 1 



2{P-x) 2 



1 duC {e \u)e^ Px = f dJ-^g(u) f dtt 2 ! dvv 2 e ivt ^ Px 

Jo s Jo Jo 

- [ duB {e \u) f dtt f dvv 2 e ivt ^ Px 
Jo Jo Jo 

(4.34) 

f duB^{u)e^ Px = f dug(u)(£ 2 + 1) f dtt 2 f v 2 dve ivt ^ Px 
Jo Jo Jo Jo 

+ 4 I duB^(u) [ dtt [ dvv 2 e lvt ^ Px 
Jo Jo Jo 

-{P-xf [ duB (e \u) [ dtt 3 [ dvv 3 e itv ^ Px 
Jo Jo Jo 



(4.35) 



i 

duA (e) (u)e l ^ Px 



o 



dug(u)(C + l) dtt 2 dvv'{i£(P-x) + l)tV 
Jo Jo 

+ 4 / duB^\u) f dtt f dvv 2 e ivt * Px 
Jo Jo Jo 



(4.36) 



Theoretically, all the chiral-even distributions B&(u), C (e) (u), A^(u) 
can be computed once the leading (physical-)twist contribution A^{u) is 
known. But, practically, it quite long and hard if one has a conformal ex- 
pansion for A^ e '(u) (except if only the asymptotic part of the wave function 
is taken, see ||). The advantage of these equations is that they are indepen- 
dent of the model used to computed (u) (another expansion can be used 
as the conformal one). They can be the basis of a numerical computation if 
one has discrete values for A^ e '(u). 



4.2 Chiral odd distributions 

For the different wave functions 

(0\u(x)a flu [x, -x]d(—x)\p~(P, A)) = 



/i / m 2 x 2 \ 

due* p " (A^(u) + ^A^(u)\ 



~\~ {P[iX v P V X ^ 



e (A) . x 



,2 r i 



1 / TYl^X^ 

//r, I due* P x ( B(°\u) + -^-B {o \u) 



o 



_L l(pW T — p( X ) T N 



P-x 



due^ p - x C (o) i 



u 



(4.37) 



the formulaes of the appendix can be used to obtains integrated equations in 
terms of the leading-twist distribution A^: 



due iiPx A {o \ 



u 



1 rl 



duE(u) dt dvv 3 e <vt * Px 
Jo Jo 



2 i£ 
1 2- 



+2 



tv 3 (Px) 3 v 3 (Px) 3 v 2 (Px) 2 
t 



v 2 (Px) 2 v(Px) v(Px) 



v(Px) 



(4.38) 
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due^ Px C (o \u) 



duS Px B {o \u 
• i 



1 /■! 



duE{u) dt dvv 2 e itv t px 
Jo Jo 



1 -4^U -2(i0 2 ^ + 2t 



_tv 2 (Px) v(Px) 



(4.39) 



(P-x) 



o 



o Jo 



1 1 z£ 

+ 



(z£)^_ 22 _t 
+ 2 1 U 2 



where 



= -z-^-A (o) (w) 
v ; 2 v ; 



tf 2 (Px) v 2 (Px) v(Px) 

(4.40) 

(4.41) 



E>( e \u) has not been computed because it is already a 0(mfy contribution. 



5 Conclusion 

This letter describes a model-independant way to extract the mass correc- 
tions of vector meson distribution amplitude from the leading twist ones. 
Although the only the (next-to-)leading twist local operators are taken into 
account, the method can be extended. In that case, one would obtain in- 
tegrated equations for meson mass corrections in function of leading and 
non-leading twist (two and three points) distribution amplitude. 

The integrated equation of this letter are quite hard to solve if conformal 
expansion model for the leading twist distribution amplitude is taken (like in 
||); but these equations can be the basis for numerical estimates, especially 
when an alternative model for the leading twist distribution amplitude is 
used (or if the distribution amplitude is directely taken form experiemntal 
results). 
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A Operators relations for different theoreti- 
cal twist 



This appendix contains some relations taken from [13], [14| and |J, which 
help to isolate the different twist parts of two-points non-local operators. 

A.l Chiral-even operator 

Consider the non-local operators ip(x)'j a ip(—x). The local expansion around 
x = can be written in the following way : 

OO j 

ip(x)-f a ip(-x) = • ..x^—i) ... D» n 7^ (A.l) 

n=0 

In order to get the leading twist contribution, the symmetrization has to 
be done over all indices and subtract the traces : 

X 

[ip(x)-f a ip(-x)] = Vi w . . .x M — ip (D w • • • -D/xn 7a) ^ - traces 

^ Ti! \ / sym 

n=0 

(A.2) 

The symmetrization has an integrated solution : 

^{x)^{-x) 

' n! I n + 1 n + 1 I 

n=0 v y 

(A.3) 

1 d - 

dv — — il)(vx)xip{—vx) (A. 4) 

OX a 

with x = 7 At x At . 

In order to have really the leading twist part of ip(x)"y a ip(— x), the leading 
twist part of ip{x)xip{—x) is needed. If one does a local expansion of this 
latter non-local operator, one will already have symetrized local operators. 
But traces has alos to be removed; this can be expressed as a differential 
equation : 

d 



dx n dx c 



x)x 



ij(-x)] lt =0 (A.5) 
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which has the formal solution : 

(x)xif)(— x)j t = i/j(x)xi/;(—x) 

oo 



oo „i 




n=l 



X' 

4 



l_f\n-l 



" (¥) 



d 2 



(n — l)!n! \dx a dx a 
Equations for the twist 3 part of ip(x) , y a ilj(— x) are the following : 



ip(tx)xip(—tx) 
(A.6) 



twist 3 



— g s / d-u / dvip(ux) uG ^ v {vx)x h '£75 



— ivG fMU (vx)x u x ip(—ux) 



+ ie„ ua P / udux v d a 



ip{x)j^ 5 ^{-x) 



twist 3 



V>(mx)7 /3 75^(-mx) 



— g s du dvip(ux) uG fJ/U (vx)x u x 



ivG fJjU (vx)x u x , y 5 ip(—ux 



+ ie,^ / udux v d a 



ip(ux)ji3il)(—ux) (A. 7) 



where is the gluon field strength, G^ = (l/2)e flua /3G al3 , and d a is the 
derivative over the total translation : 



ip{ux)^pip{—ux) 



_d_ 

dy a 



ip{ux + y)jpip(— ux + y) 



(A.8) 



y-*0 



These equations come from [[14]] and can be obtained with the equation 
of motion, neglecting quark masses. 

Another useful formula is the equation (5.13) in which can be written 
like this : 







dx n dx a 



ip(x)xip(—x) = —d 2 ip{x)xip{—x) + operators with gluons (A. 9) 



The expansion of [A.6| at the order x , gives 



ip(x)xip(—x) 



x)x^(-x)] l t + -x 



4 
12 



dv- 



d 2 - 



dx n dx 



—ip{vx)xip{—vx) + . . . 



x)xifj(—x)\ l t x / dvv d i/j(vx)xip(—vx) 



+ operators with gluons + twist 3 + . . . 



(A.10) 



This last equation gives a twist 4 contribution with a total derivative 
(defined in |A.8| ). 



A. 2 Chiral-odd operator 

Consider the non-local operators ^(x)a a /3ip(— x). The following equation can 



be used[17 



i/j(x)a a pip(-x) 



1 d - 

udv — — ip(vx)a au x u ip(~ vx ) — (a <-> /3) + 0(twist 3) 



o 



(A.11) 



to extract the most symetrical part of ip(x)a a p^(— x): 



ip(x)a al3 '4)(-x) 



J sym 



1 d - 

vdv — — ip(vx)a au x u ip{—vx) — (a «-> /3) 



dx 



P 



(A.12) 



(it reproduces the result in [plj] ). 

In 0, the equations for the leading twist part of %{){x)a au x l 'ip{—x) can be 
founds : 



d 
dx a 
d 2 
dx p dx p 



x)a ap x 1J i)(-x)] x ^ = 
x)a av x v i){-x)] L t _ = 



(A.13) 
(A.14) 



Like in equation A. 10, the contribution of total derivative can be added : 



ip{x)a au x v ip{—x) 



x)a au x u ip{~ x)] j t — -a; 2 ^^ / dtt 2 d 2 (ip{tx)a au il ! {—tx)) 



+ operators with gluons + . . . 



(A.15) 
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